Quantum stabilization of 1/3-magnetization plateau in Cs 2 CuBr 4 
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We consider the phase diagram of a spatially anisotropic 2D triangular antiferromagnet in a magnetic field. 
Classically, the ground state is umbrella-like for all fields, but we show that the quantum phase diagram is much 
richer and contains a 1/3 magnetization plateau, two commensurate planar states, two incommensurate chiral 
umbrella phases, and, possibly, a planar state separating the two chiral phases. Our analysis sheds light on 
several recent experimental findings for the spin-1/2 system Cs2CuBi'4. 

PACS numbers: 



Introduction. A defining characteristic of frustrated quan- 
tum magnets is the appearance of numerous competing or- 
ders. This competition dramatically enhances quantum fluc- 
tuations, generating highly non-classical behavior as exem- 
plified by, e.g., CS2CUCI4 and CS2CUB14. These materi- 
als comprise quasi-2D spin-1/2 triangular antiferromagnets 
with spatially anisotropic exchange [see Fig. [TJ a )] an d weak 
Dzyaloshinskii-Moriya (DM) coupling. Absent the latter, 
both systems classically should realize a zero-field coplanar 
spiral, which evolves into non-coplanar "umbrella" states in a 
field as in Fig. [TJb) with smoothly increasing magnetization 
up to saturation [1]. Experiments, however, reveal decidedly 
different, non-classical behavior: in fields directed along the 
triangular layers CS2CUCI4 realizes commensurate coplanar 
order in a wide field range with smoothly increasing magneti- 
zation 12j,|3[], and Cs2CuBr4 exhibits collinear "up-up-down" 
(UUD) order shown in Fig. [TJc) over a finite field interval, 
yielding a 1/3 magnetization plateau Bl S 01- Neither 
observation is accounted for within a classical analysis JH. 
NMR 1 8, 9] and neutron scattering [6] additionally find planar 
states adjacent to the UUD phase, with neutron and thermo- 
dynamic measurements J3l indicating that the transitions are 
first order. Additional experiments fla il (XI on Cs2CuBr4 also 
suggest the presence of a narrow 2/3-plateau and additional 
intervening collinear phases near particular fields. 

While the existence of the UUD phase is well-established 
for the isotropic triangular antiferromagnet, much less is 
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FIG. 1: (a) Anisotropic triangular lattice with horizontal exchange 
J and diagonal exchange J', (b) Umbrella and (c) planar phases 
comprise competing classical ground states of the isotropic nearest- 
neighbor model. 
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FIG. 2: Proposed phase diagram for the anisotropic nearest- 
neighbor Heisenberg model near 1/3 magnetization (full field range 
not shown). The horizontal axis is 8 = (40/3)5(7 - J') 2 /J 2 . Pla- 
nar states shown are commensurate, though they are expected to be 
incommensurate at small and large fields. The shaded area is where 
the UUD and adjacent phases are metastable, the energy being mini- 
mized by umbrella states of Fig. [TJb). 



known about the stability of the plateau and the proximate 
quantum phases in the anisotropic case. The challenge here 
is illuminated by first observing that in the isotropic limit, the 
UUD state appears due to an "accidental" classical degener- 
acy between umbrella and planar states shown in Figs. [TJb) 
and (c), which quantum fluctuations lift in favor of the latter 
lilUl . When J' 7^ J, however, this degeneracy is lifted already 
at the classical level, but in favor of umbrella states for all 
fields. The planar phases can then only emerge if quantum ef- 
fects overshadow those of spatial anisotropy. This in turn im- 
plies that the standard spin-wave expansion is not applicable 
since planar phases cease to be classical ground states. To ad- 
dress the quantum phase diagram for the anisotropic system, 
particularly the phases near 1/3 magnetization, we introduce a 
modified approach here which is controlled by the smallness 
of 1/5 and spatial anisotropy, and yields results which are 
non-analytic in both parameters. 

Figure[2]summarizes our results. We find that the physics is 
controlled by the parameter S = (40/3)5(7- J') 2 / J 2 . For 
S < 1 the stability of the UUD phase is, counter-intuitively, 
unaffected by anisotropy. Moreover, the spin order remains 
coplanar and commensurate at fields both below and above 
the UUD phase; incommensurate phases appear only at small 
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and high fields. For 1 < 6 < 4, the UUD phase per- 
sists, but at the boundaries it becomes unstable towards non- 
coplanar, incommensurate phases which can be regarded as 
distorted umbrellas (this happens for 6 > 1 at the lower 
boundary and for 5 > 3 at the upper boundary). These 
two phases emerge as finite-fc instabilities of the two low- 
energy spin-wave branches of the UUD phase, and both have 
a non-zero Ising order parameter associated with chirality 
Kabc = z ■ {S A x S B + S B x S c + S c x S^) for each pla- 
quette. For 6 > 4, the UUD state ceases to exist, and there is 
no magnetization plateau. Since the chiralities of the low- and 
high-field distorted umbrella phases are uncorrelated, the two 
must be separated by (at least) a first order transition in this 
region. At still larger S (stronger anisotropy), ID physics be- 
comes important, and the system cannot be described by our 
semi-classical theory. 

As a further complication, for 6 > 2 the energy of the UUD 
state becomes larger than that of the classical, undistorted um- 
brella, i.e., for 2 < 6 < 4, the UUD state and neighboring dis- 
torted umbrellas are metastable. We represent this by shading 
the region 6 > 2 in Fig. We expect that these metastable 
phases may be probed in pulsed field experiments 11211 . 

Model and UUD state in the Anisotropic System. We con- 
sider a simple Heisenberg model with 
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where S r are spin-5 operators, the exchanges J rr > are as 
shown in Fig. Ufa), and h is the (scaled) magnetic field. The 
saturation field is h sa ,t = (2J + J') 2 / J. Since we wish to 
treat quantum effects and the effects of anisotropy on equal 
footing, we will organize our analysis by assuming that both 
(J - J')/ J and 1/5 are small. 

With J = J', the two competing classically degenerate 
states are commensurate (three-sublattice) umbrella and pla- 
nar states shown in Fig. Q] Quantum fluctuations favor pla- 
narity, and spin re-arrangement in a field occurs as in Fig.QJc). 
This process includes an intermediate UUD phase, which is 
classically stable only at h sat /3, but quantum fluctuations ex- 
tend its stability to a finite field interval, hP cl < h < hP c2 Jllll . 
resulting in a 1/3-magnetization plateau. This is not surpris- 
ing given that quantum fluctuations generally favor collinear 
states {HQ. To leading order in 1/5 
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which for 5 = 1/2 yields a plateau in a range Ah — 
h° c2 — hP cl = 1.8 J I (25*), in good agreement with exact diag- 
onalization [15]. Inside this range, there are two low-energy 
spin-wave modes with gaps oc \h° cl 2 — h\ at k = 0. 

When J J', the umbrella state becomes incommensu- 
rate, and classically has lower energy than the planar phase 
for all fields. The naive expectation, then, is that the UUD 
phase must immediately shrink and disappear as | J — J'\ in- 
creases. We show, however, that the actual situation is much 
more complex, with new phases emerging when J ^ J'. 



To study the stability of the classically unfavorable UUD 
state, we explore a modified large-5 approach to Eq. ([T). First, 
we introduce a three-sublattice representation where spins on 
the A and B sublattices point up while those on the C sub- 
lattice point down, and use the standard Holstein-Primakoff 
mapping. The usual linear spin- wave Hamiltonian obtained in 
this fashion is not an appropriate starting point due to the clas- 
sical instability of harmonic spin waves at 6 ^ 0. However, 
the interacting spin-wave Hamiltonian must support a stable 
UUD plateau over a finite anisotropy range, as exact diagonal- 
ization finds II 1 811 . Therefore, we extend the linear spin-wave 
Hamiltonian of the UUD state to include the leading 1/S self- 
energy corrections obtained by decoupling the quartic interac- 
tions using correlations from the isotropic system. 

Diagonalizing this Hamiltonian, we obtain three spin-wave 
branches. One branch describes a precession of the total mag- 
netization and has a high energy ~ ho = J + 2 J'. The ener- 
gies of the other two branches are small near k = 0: 

-Huud = 5^[cjidj k di, k + w 2 4.k rf 2,k], (3) 

k 

where the leading expressions at small k are 

u**0_ ± (»-».--^ -?*»)+ 1§ (4) 



with Z = ^9 + 105[6fe2 + 105fc4 - 3k*{6 - 2)]. The criti- 
cal fields obtained from these energies are 



3 J / Z — 3 

h c i, 2 = h° cl2 + 2(J' - J) T — min l T k£ + 
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where h° cl , 2 is given by (0, and the minimum is taken with 
respect to k x (k y =0 at the minimum for all S). The UUD 
phase is stable for h c \ < h < h c2 . 

These results, which are non-analytic in 1/S and J — J', 
encode the physics governing the local stability of the UUD 
state in the anisotropic system. One can verify by sending 
S — > oo above that the UUD state is indeed unstable for any 
non-zero anisotropy in the classical limit, due to an instabil- 
ity at finite k x . Surprisingly, in the quantum system a finite 
amount of anisotropy is required to begin destabilizing the 
plateau. Specifically, for 5 < 1 both modes are minimized 
at k — 0, so it follows from Eq. (f5]) that the plateau width 
Ah is unchanged from the isotropic system. The effect of 
anisotropy in this regime is only to shift the plateau's location 
and soften the dispersion around k = 0. 

For S > 1, the minimum of oji shifts to ki± = (±£;i,0), 
where k\ = [36 - 6 + y/3S(A - 6)]/(20S); the lower criti- 
cal field then moves upward, reducing the width of the UUD 
plateau (see Fig. |2). Similarly, for 5 > 3 the minimum 
of uj 2 shifts to k 2 ± = (±fc 2 ,0), with k 2 = [35 - 6 - 
a/35(4 - <J)]/(205). At this point the upper critical field 
moves to a smaller value, further reducing the UUD region. 
The plateau ceases to be locally stable at S = 4, when both 
spin-waves become gapless at k\ = k\ = kf n = 3/(105). 
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Let us now explore the phases that emerge immediately 
away from the UUD state. At h c \ and h c i, magnons Bose 
condense, and one must determine the energetically favorable 
combination of operators <ix,2k tnat condenses, and what this 
implies for the spin components (S x,v ). For 8 < 1 this is 
straightforward: the minima of wi^k) occur at k = , and 
the order parameters are simply ^x,2 ^ (di,2,o)- One can eas- 
ily verify that condensation of ipi(ip2) at h — h c i (/i C 2) leads 
to the commensurate coplanar spin configurations displayed 
in Fig. Q2 C )- The prediction of commensurate order adjacent 
to the UUD state over a range of anisotropy is rather nontriv- 
ial, and could be tested in exact diagonalization studies. 

The situation is subtler at the lower critical field when 
8 > 1, since here Wx(k) possesses two inequivalent min- 
ima at ki±. There are then two order parameters, ip± = 
y/3/NS(di^ 1± ) {N is the number of spins), whose energy 
derived from the interacting spin wave Hamiltonian IU6I1 is 
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Here r oc h — h c \ and u = 2 cosh 2(j>k 1 > where 
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Since u > 0, below the transition interactions fa- 
vor tp + 0, ip- — or vice versa. Choos- 
ing the former, the spin configuration can be writ- 
ten (S\) = — iS^+fcosh^fej + ismh<j> kl )e~ lklX , 
(S+) = Sty+(cosh0 fcl + isinh0 fel )e- lfc ^, (S+) = 
2iS4>+ smh(f>k 1 e + ' lklX . This corresponds to non-coplanar, 
incommensurate order that can be described as a distorted 
umbrella. Non-coplanarity of this state leads to a finite 
chirality K^ l \ the sign of which is determined by that of the 
condensate momentum via K% c = ±3S' 2 |V'±| 2 sinh20 fcl . 

The same consideration holds at the upper critical field 
when 8 > 3: ^(k) again has two inequivalent minima at 
k2±, and the energy has the same form as in (O, with the 
order parameter ip± = \/3/NS(d2,k 2 ±)- The spin configu- 
ration above h C 2 is another distorted umbrella with chirality 
^i 2 ic = T3S 2 |Vi ± | 2 sinh20 fe2 . 

At 8 = 4, the UUD plateau shrinks to a point at h c = 
ho + 17 J /(AOS), and becomes unstable at larger 8. How the 
two distorted umbrellas merge in this regime presents an in- 
teresting issue. Since these states arise upon condensation of 
different spin-wave modes at 7i c i,2. their chiralities are uncor- 
rected. The two phases then cannot gradually transform into 
each other and must be separated either by a first order transi- 
tion, or by an intermediate phase with no chirality. 

To gain insight here we study the instability of the UUD 
phase at 8 = 4, h = h c . At this point, the two spin-wave 
branches become gapless at the same ±fc m , and the coher- 
ence factors smh(f>k m and cosh</>fc m diverge as 1 /V4 — 8, so 
that tanh2(/>fc m — > 1. There are more choices for the order 
parameter ath = h c compared to either h c \ or h C 2 as both ip± 
and tjj± condense at h c . The full expression for the ground 



state energy at 8 = 4, h = h c to fourth order in ip and ip, and 
to leading order in 1/S is 
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subject to constraint ip—ipl. — ip^ip + — i(ip—ip_ + V'+' ! /'+) 
which eliminates infinitely large terms from the energy. 
Choosing just one of the four order parameters non-zero, we 
obtain the same distorted umbrella states as before, with E oc 

4 and a finite chirality. However, we see from ([8]) that there 
is a better choice — taking = 7^ 0, ip- = t/j + = 
or vise versa, we find that E = 0. Thus, unlike the situation 
at any other point on the critical lines h c i,2(8), the magnitude 
of the condensate = at the end-point of the UUD 
phase is unconstrained implying that |^ + | = |-0_| jumps to 
a finite value right at 8 = 4. The chirality Kabc of such a 
state depends on the relative phase 8 of the two order param- 
eters as Kabc oc tanh 2(j>k m + sin 9 and vanishes when 8 — 
— arcsin(tanh 2<pk m ) — * — t/2. We verified that this particu- 
lar 9 is the only choice at which the transverse magnetization 
given by (S x /) - 0, (S+) = -(S+) ~ \^+\V^S e~ ik ^ 
does not diverge together with the coherence factors but rather 
remains zero at 8 = 4, even though ^ there. It is 
tempting to speculate that such a zero-chirality state persists 
beyond 8 — 4 along a line in the h — 8 plane and continues 
to separate the two distorted umbrella phases (it cannot ex- 
ist in a finite /i-range since, unlike the UUD phase, it does 
not have two gapped low-energy modes). The spin structure 
along this line is either collinear, as at 8 — 4, or coplanar, 
with (S X X V ) = 0, (5£) = -(Sc); the difference can not be 
resolved within our formalism. 

Energy considerations and the phase diagram. So far we 
have analyzed the UUD phase's local stability without ad- 
dressing whether it globally minimizes the energy. There are 
three regimes where one can easily compare the umbrella and 
planar energies. First is the high field regime h w /i sat . There, 
the umbrella state, which at arbitrary h is described by 

S r = ^{cos 6>[cos(Q • r)x + sin(Q • r)y] + sin 8z} (9) 

with Q = 2 cos _1 (— J'/2 J) and sin# = h/h sat , wins for 
all J' 7^ J simply because quantum effects vanish at h sat . 
We have verified this explicitly by computing the analog of 
Eq. © at the saturation field to show that indeed interactions 
drive the system into the umbrella state for arbitrary J' ^ J. 
As a result, the critical line which begins at 8 = 3, h = h C 2 
should end up at 8 = 0, h = h sat . 

The second regime occurs at small h — > 0. Here the lowest- 
energy planar configuration is incommensurate, with the same 
Q as the umbrella state and 
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sin(Q ■ 



(10) 



where <p T — —(2h/u) sin(Q ■ r) + 0(h 2 ) and u = h sat [l + 
(J — J') 2 1 J 2 ]. At small h, the energy difference between the 
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incommensurate umbrella and planar states of Eqs. (0 and 
m is AE h ^ = (E umh - E pl )/NS 2 = -(l/2)/i 2 A x , 
where Ax = Xumb — Xpi is the difference of susceptibili- 
ties. In the classical limit we find Xumb = l/^sat) Xvl = l/ u > 
so that A% = (J - J') 2 /(9J 3 ) and the umbrella state has 
lower energy. The competition comes from quantum fluc- 
tuations: 1/5 corrections to Xumb an d Xpi are different al- 
ready for J = J', and such that Ax qu ~ -0.16/(18J5) 
(Ref. [11]). Adding the two contributions, we find that 
AE h ^„ = [0.008/i 2 /(2J5)](l.l - 6), i.e., the incommen- 
surate planar state has lower energy for S < 1.1. This implies 
that the commensurate planar state that we found immediately 
below hd should undergo either a second- or first-order tran- 
sition into an incommensurate planar state at some h < h c i- 
We therefore expect the line separating planar and distorted 
umbrella states at low fields to depart at 6 = 1, h = h c \ and 
end up at S = 1.1, h = 0. 

Finally, at /i sa t/3 the energy difference between the um- 
brella and UUD phase is AE 1/3 = [0.067 J/(25)] (2.0 - 6), 
where the first and second terms, respectively, are the classi- 
cal and quantum contributions; see Ref. [11]. Consequently, 
the UUD phase and the neighboring distorted umbrella phases 
remain global minima only up to 5 = 2.0 and become 
metastable at larger 5. This suggests that for 5 > 2 the 
UUD state can be observed only via a transient magnetiza- 
tion plateau, similar to the situation in a kagome system 11211 . 
Equilibrium measurements should reveal only umbrella-like 
states in that region of 6. 

The resulting phase diagram near 1/3 magnetization is 
shown in Fig. [2] It contains an UUD phase; two commen- 
surate planar states from Fig. [TJ C ); an d two non-coplanar in- 
commensurate distorted umbrella phases. The shaded region 
corresponds to the regime where the classical umbrella mini- 
mizes the energy globally. Additionally, incommensurate pla- 
nar states are expected at small fields when S < 1, and near 
the saturation field for small 6. We also expect new phases 
at small J'/ J (large 5), where one-dimensional physics takes 
over the semi-classical analysis. 

This phase diagram is in agreement with data for 
Cs2CuBr4, where J'/ J = 0.7 implies that 6 = 0.6 if we 
extrapolate to 5 = 1/2. For this 5, the UUD state is present, 
and the nearby phases are planar, in agreement with NMR 
|JsL 13] and neutron |(J experiments. These experiments also 
observe that both transitions out of the UUD state are first or- 
der. Our calculations predict continuous transitions as a con- 
sequence of the U(l) spin symmetry exhibited by the Hamil- 
tonian (HJ. However, when this U(l) symmetry is broken 
explicitly by spin-orbit coupling, cubic terms in the free en- 
ergy are permissible, which generically render the transition 
first order. In particular, DM coupling of the form present 
in Cs2CuBr4 breaks this symmetry when the field is directed 
along the triangular layers. In addition, a direct first order 
transition from UUD phase into the incommensurate planar 
phase is also a possibility, which should be investigated by nu- 
merical calculations similar to those in 11711 . For CS2CUCI4, 
the anisotropy is much higher (S « 2.9), and the system very 



likely lies outside of the applicability region of our analysis, 
and should be approached from a ID perspective [19fl . Still, 
even within our framework, 8 > 2 implies no UUD phase, 
and no plateau is seen in CS2CUCI4. 

An intriguing question concerns the possible appearance of 
a 2/3-magnetization plateau at <5 < 1 and higher fields, as ob- 
served in Cs2CuBr4[|5[], which would correspond to (at least) 
a "5-up, 1-down" configuration. While such states are never 
ground states to order 1/5, we verified that their energy is 
reduced when J' 7^ J. We speculate that, due to a large de- 
generacy of 5-up, 1-down configurations, a 2/3-plateau may 
be entropically stabilized at finite temperature. Regarding this 
issue, the role of spin-phonon couplings should be seriously 
investigated 12011 . 

Conclusions. Using a modified large-5 approach, we stud- 
ied the quantum phase diagram of an anisotropic triangular 
antiferromagnet, with particular emphasis on the classically 
unstable UUD state and proximate phases. Fig. [2] summarizes 
our findings. The UUD phase with 1 /3 magnetization plateau 
survives a substantial range of anisotropy, and at its bound- 
aries transforms either into commensurate planar phases, or 
into umbrella-like incommensurate chiral phases, depending 
on whether the spin-wave instabilities of the UUD phase are 
at zero or finite momenta. Our results explain a number of 
experimental findings for Cs2CuBr4. 
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